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1. Introduction 



Abstract. Using the framework of the algebraic Bethe Ansatz, we study the scalar 
product of the inhomogeneous XXZ spin-^ chain. Inspired by the Izergin-Korepin pro- 
cedure for evaluating the domain wall partition function, we obtain a set of conditions 
which uniquely determine the scalar product. Assuming the Bethe equations for one set 
of variables within the scalar product, these conditions may be solved to produce a de- 
terminant expression originally found by Slavnov. We also consider the inhomogeneous 
XX spin-i chain in an external magnetic field. Repeating our earlier procedure, we find 
a set of conditions on the scalar product of this model and solve them in the presence of 
C " 3 ■ the Bethe equations. The expression obtained is in factorized form. 

04 , Scalar products are objects of essential interest in models solvable by the quantum inverse 

scattering method/algebraic Bethe Ansatz [T]. In the case of two-dimensional integrable 
models, the definition of the scalar product may be phrased as follows: Let Va,Vb denote 
' copies of the vector space C^, and consider a generic two-dimensional model with i?-matrix 

Oh, Rab{u,v) G End(Va Vb), monodromy matrix 

and pseudo- vacuum states (0|, |0) satisfying 

(2) {0\A{u) = a(M)(0|, {0\B{u) = 0, (0|£'(m) = diu){0\, 

(3) A{u)\Q) = a(u)|0), C(u)|0) = 0, D{u)\0) = d{u)\0), 
m ; (4) (0|0) = 1. 

The scalar product of the model is then given by 

(5) Sn[{u},{v}) = {0\C{ui) . . .C{un)B{vi) . . . B{vNm, 



(N 
^' 

■ where the operators C{u) and B(v) obey a non-trivial commutation relation which is recov- 

ered as a single component of the bilinear equation 

(6) Rab{u, v)Ta{u)Tb{v) = Tb{v)Ta{u)Rab{u, v) . 

_ Finding simple and exact expressions for scalar products is important for at least two reasons: 

rS ■ 1- In the limit when the variables {u} and {v} are equal to the same solution of the Bethe 

I equations, they reduce to the norm-squared of a Bethe eigenstate; 2. Algebraically, they 

" ■ ■ are closely related to correlation functions and constitute the first step in the calculation 

of such quantities. Nevertheless, the explicit calculation of the scalar product is a difficult 
task, owing to the complicated nature of the algebraic relations contained in 1^ . 

A major development in the theory came in [2], when N A Slavnov found a determinant 
expression for the scalar product ([S]) of models with an A^^^^ type i?-matrix [3] in the case 
when one set of variables {v} satisfies the Bethe equations, whilst the other set {u} is left 
free. Aside from its intrinsic beauty, the Slavnov scalar product formula has many practical 
uses. The simplest of these is seen by setting the free set of variables {u} equal to the set 
{u}, when it specializes to a formula for the norm-squared of a Bethe eigenstate, originally 
proposed in [i] and proved in [5]- The Slavnov formula was also indispensable to the work 
of [6], in which a compact expression was found for a generating functional of correlation 
functions in generic A^^-* type models. More recently, by virtue of its determinant form, it 
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was shown that the Slavnov scalar product can be viewed as a specialization of a r-function 
of the KP hierarchy ^ or more simply as a solution of discrete Hirota equations j8]. 

The main purpose of this paper is to provide a recursive, diagrammatic proof of the 
Slavnov formula in the particular case of the inhomogeneous XXZ spin-^ model. In order 
to do this, we resort to the relationship between the XXZ model and the six-vertex model 
of statistical mechanics. We write the XXZ scalar product as a sum over lattice states in 
the six-vertex model, subject to certain boundary conditions. The lattice sum is shown to 
satisfy a set of four conditions which determine it uniquely, and then the Slavnov formula 
is proposed as a solution. This type of procedure is motivated by the work of V E Korepin 
and A G Izergin on the domain wall partition function of the six-vertex model, whereby a 
set of conditions was found in [S] and subsequently solved in [S]. 

In [TU], B U Felderhof defined a transfer matrix which commutes with the Hamiltonian 
of the XY spin-i chain in an external magnetic field. The transfer matrix thus constructed 
was parametrized in terms of Jacobi elliptic functions. By specializing to the trigonometric 
limit, one obtains a transfer matrix which commutes with the Hamiltonian of the XX spin-^ 
chain in an external field. We shall consider this latter case in the second part of the paper, 
and refer to it as the trigonometric Felderhof model. It previously appeared in [TT] as the 
first in a hierarchy of vertex models, and its domain wall partition function was calculated 
in [2] using an Izergin-Korepin type of procedure. In reflection of the free-fermionic nature 
of this model the domain wall partition function has a product form, in which all zeros 
clearly factorize. In this work we calculate the scalar product of the trigonometric Felderhof 
model, once again assuming the Bethe equations for one set of variables. We find that the 
scalar product also has a product form, with all zeros factorized. Some of the zeros are of 
the type found in the domain wall partition function, while the rest are roots of the Bethe 
equations. 

The outline of the paper is as follows: In §2 we introduce basic facts related to the 
XXZ spin-i chain, describing its solution via the algebraic Bethe Ansatz. In §3 we define 
the domain wall partition function of the six-vertex model, and review the method of its 
evaluation due to [SUS]. This serves as a useful introduction to §4, where we present our 
Izergin-Korepin procedure for calculating the scalar product. We set up a sequence of 
intermediate scalar products {S„}o<„<n which were originally considered in [13], and write 
down conditions which uniquely determine them, before proposing a solution and showing 
that it is valid. In §5 we introduce the XX spin-^ chain in an external magnetic field, and 
describe the method of diagonalizing its Hamiltonian due to Felderhof. The calculation 
of the trigonometric Felderhof domain wall partition function is reviewed in §6, before we 
evaluate its scalar product in §7. Once again, an Izergin-Korepin approach is employed in 
calculating these quantities. 



2. XXZ SPIN-i CHAIN 

In this section we introduce the basics of the XXZ spin-i chain, describing its solution 
via the quantum inverse scattering method/algebraic Bethe Ansatz. For a more general 
introduction, the reader is referred to chapters VI and VII of the book [T]. 

2.1. Space of states V. The finite length XXZ spin-i chain consists of a one-dimensional 
lattice with M sites. Each site m is occupied by a spin-^ fermion and has a corresponding 
two-dimensional vector space Vm with the basis 

(7) Basis(V;„) = { tm,;™ }, 

where for convenience we have adopted the notations 

Physically speaking, t„i and represent the spin eigenstates of a spin-^ fermion at site m. 
The space of states V is defined as 

(9) V = Vi<E)---<»Vm 
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and it is the goal of the algebraic Bethe Ansatz to find states within V which are eigenvectors 
of the Hamiltonian (see . The two simplest states in V are those for which all spins are 
up/down, and we prescribe them the notation 

M M 

(10) I i\m) = (g) u, I ^m) = (g) ;™ . 

rn—l m—1 

We will also make use of the definitions 

(11) \1\N/Al)^ (g) tm (g) U, I^AT/Af)- (g) U (g) tm 

l<m<Ar N<m<M l<7n<N N<m<M 

for states whose first N spins are up/down, with all the remaining spins being down/up, 
respectively. 

Analogous definitions apply in the construction of dual spaces. To each site m we associate 
the dual vector space with the basis 

(12) Basis(y,;) = {t:„,;:„}, 

where we have adopted the notations 

(13) t:„ = ( 1 ),„, = ( 1 

and from this we construct the dual space of states 

(14) V* ^ V^* (g> ■ ■ ■ (g> V^f. 
Similarly to before we shall write 

M M 

(15) (^M| = (gt:„, {Ul\ = (^i*m 

m—1 m—1 

for the dual states whose spins are all up/down, and 

(16) {i[N/M\= (g t:„ (g (^at/mI^ (g ;™ (g t™ 

l<m<Ar N<m<M l<m<N N<m<M 

for dual states whose first N spins are up/down, with all remaining spins being down/up, 
respectively. 



2.2. Pauli matrices. We define the Pauli matrices 
(17) < - 



1 \ ,y ^ I 

10/'" \ i / ' " I -1 



with i — and the spin raising/lowering matrices 

(18) = + = ) ' " h^"""' ^ ''^'"^ " ( 1 
where in all cases the subscript m is used to indicate that the matrices act in the space Vm- 

2.3. Hamiltonian H. The Hamiltonian of the finite length XXZ spin-i chain is given by 

(19) ^ = ? E + «+i + A(«+i + 1) 



m—1 



where A = cosh 7 is the anisotropy parameter of the model, and the periodicity conditions 
"'m+i = '^ii'^'m+i ~ '''i''''if+i = ^1 ^'^^ assumed. In the coming subsections, we shall 
review the algebraic Bethe Ansatz procedure for finding the eigenvectors \^) d V of H . 
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2.4. i?-matrix, crossing symmetry and Yang-Baxter equation. The i?-niatrix corre- 
sponding to the XXZ spin-i chain is given by 



(20) 



Rab{u,v) 



/[u-u + 7] \ 

[u~v] [7] 
[7] [u-v] 

V 00 [u-v + j]J^^ 



where we have defined [u] = sinhu. The i?-matrix is an element of End(Va ® Vf,), where 
Va, Vb are copies of C^. The variables u,v are rapidities associated to the respective vector 
spaces Va,Vb, and 7 is the crossing parameter. We identify the entries of (00)) with vertices, 
as shown in figure [T] 




{u^ v) — [u — V + 7] 



fc+ {u, v) — [u — v] 




a_ (u, d) — [ti — i; + 7] 



b-(u, v) — [u — v] 




Figure 1 . Six vertices associated to the XXZ i?-matrix. Each entry of the 
i?-matrix (l20l) is matched with a vertex. All types of vertex not shown are 
by definition weighted to zero. 



Lemma 1. Define u = it + 7 for all rapidities u. The i?-matrix has the crossing symmetry 
property 



(21) 



Rab{u,v) = ~(TlRba[v,uf''(jl 



where cr^ is the second of the Pauli matrices (I17p acting in VJ,, and t;, denotes transposition 
in the space End(Vf,). 



Proof. We have 

(22) Rba{v,uY^^ 
which leads to the equation 

(23) -alRba{v,uf^al 



( [v — u-\- ^] 



V [7] 



[7] \ 

v-u\ 

[v-u] 

[v-u + -f] ) 



ha 



I 



[u — w — 7] [7] 



\ 







V 



[7] [u — t; — 7] 





Finally, using the definition u = u + 7 and the fact that Rba{u,v) = Rab{u,v), we prove 



□ 



Lemma 2. The i?-matrix obeys the Yang- Baxter equation 

(24) Rab{u, v)Rac{u, w)Rbc{v, w) = Rbc{v , w)Rac{u, w)Rab{u, w), 

which holds in End(Va Vf, (g) K) for all m, i;, ui. 



2.5. L-matrix and local intertwining equation. The i-matrix for the XXZ model de- 
pends on a single indeterminate u, and acts in the space Va- Its entries are operators acting 
at the m}^ lattice site, and identically everywhere else. It has the form 

(25) L^,n{u)=( [7^^™] . 
where we have defined, as before, [u] — sinhu. 

Lemma 3. Using the definition of the i?- matrix ^U\i and the L-matrix (pS)) . the local 
intertwining equation is given by 

(26) Rab{u, v)Larniu)Lb„i{v) = Lbmiv) Larn{u) Rabiu, v). 

Proof. Using the matrix representations of {c™, cr" , cr^} given by equations (jl7p and (jlSp . 
we find that the L-matrix (1^51) takes the form 



(27) Lam(u) 



/ [u+"j] 

[u-i] [7] 

[7] [^^-i] 

V [u+^] J 



-Ram (U, 7/2), 



from which we see that it is equal to the i?-matrix Ram{u,w„i) with Wm = 2- local 
intertwining equation (j26p becomes 



(28) Rab{u, v)Ram{u, 7/2)i?b„(w, 7/2) = Rb„i{v, 7/2)i?am(w, j/2)Rab{u, v), 

which is simply a corollary of the Yang-Baxter equation ([M]). 



□ 



2.6. Monodromy matrix and global intertwining equation. The monodromy matrix 
is an ordered product of L-matrices, taken across all sites of the chain. It is given by 

(29) Ta{u) = Lal{u)...LaM{u). 

In the last subsection we observed that the L-matrix Lam{u) is equal to the i?-matrix 
Ram{u,Wm) uudcr the specialization Wm = -^^ Using this observation, it is convenient to 
construct an inhomogeneous monodromy matrix as an ordered product of the i?-matrices 
Ram{u,Wm)^ without restricting the variables Wm- That is, we define 

(30) Ta{u, {w}m) = Ral{u,Wi) . . . RaM {u,W m) ■ 

The variables {wi, . . . , wm} are called inhomogeneities and the usual monodromy matrix is 
recovered by setting Wm = ^ for all 1 < m < M. It turns out that the inclusion of the 
variables {wi, . . . ,wm} simplifies many later calculations. The contribution from the space 
End(T4) can be exhibited explicitly by defining 

/oiA T S,n\ \-( Mu,{w}m) B{u,{w}m)\ 

V C'(m,{w}m) D{u,{w}m) ) ^ 



where the matrix entries are all operators acting in V 
version of these operators is given in figure [21 
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® Vm ■ The diagrammatic 



A{u, {w}m) 



S(u, {w}m) 



-t> 



IDl Wm 

C{u,{w}m) 



wi Wm 

D{u, {w}m) 



Figure 2. Four vertex-strings of the XXZ monodromy matrix. Each entry 
of (|3ip is matched with a string of i?-matrix vertices. 

Owing to the Yang-Baxter equation (|24p and the definition of the monodromy matrix ([SO)), 
we obtain the global intertwining equation 

(32) Rab{u,v)Ta{u,{w}M)Tb{v,{w}M) = Tb{v,{w}M)Ta{u,{w}M)Rab{u,v) 

This equation contains sixteen commutation relations between the entries of the monodromy 
matrix. Particular examples of these commutation relations include 

(33) B{u)B{v) = B{v)B{u), 

(34) [u~v + j]B{u)A{v) = [j]B{v)A{u) + [u - v]A{v)B{u), 

(35) [-f]B{u)D{v) + [u- v]D{u)B{v) ^ [u - v + -f]B{v)D{u), 

(36) C{u)C(v) = C(v)C(u), 

(37) [-f]A{u)C{v) + [u- v]C{u)A{v) ^[u-v + j]A{v)C{u), 

(38) [u~v + j]D{u)C{v) = [j]D{v)C{u) + [u ~ v]C{v)D{u), 

(39) D{u)D{v) = D{v)D{u), 

which are obtained by multiplying the matrices in (1321) and equating components. 

Lemma 4. Let {u} l — {ui, . . . , u^} and {w} m — {wi, . . . , wm} be two sets of variables, 
with cardinalities L,M > 1. To each variable ui we associate an auxiliary space VJj, , while 
to each variable Wm we associate a quantum space Vm- Defining 

(40) t(^{u}l,{w}m^ ^ Ta^{uL, {w}m) ■ ■ - Ta^iui, {w}m), 
we claim that 

(41) {w}m) = (-)^*'Ti(w;i, {u}l) . . .Tm^wm , 
where for all 1 < m < Af we have defined 
(42) 

with {u}l ^ {ui +J, . . . ,ul + 7}. 
Proof. By the definition (HUl) we have 
(43) 

t({u}l,{w}m^ = (^RaLl{uL,Wi) . ■ ■RaLM{uL,WM)j ■ ■ ■ (^ail(ui, Wl) • . . i?ai M ("l , WM )) • 

Commuting i?-matrices which act in different spaces leads to the equation 

(44) 

t({u}l,{w}m^ = (^RaLliuL,Wi) . . .i?ail(wi, Wl)) • • • (^Ra^ M {ul, Wm) ■ ■ ■ RaiM{ui,WM)^ , 
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m{Wm, Wl) - I _C{w^,{u]l) A{Wra. {u} l) 



and using the crossing symmetry relation (j2ip on every i?-matrix in (j44p we obtain 
(45) t(^{u}lAw}m) <jl{Riad^i,ULT' ■■■RiaAwuUiY')al 

X ...X alj {RMaL {WM-, UlY" . . . RMai {wm-, 



From the standard identity of matrix transposition {AB)* — B^A* we may reverse the order 
of the i?-matrices in (l45t . yielding 

(46) t({u}l,{w}m) = (-)^*Vf(i?iai(l«l,Sl)...i?laJu'l,Si))*V^ 

/ _ _ \ tAf 

X ...X (T%i iRMai iwM,Ui) . . . RMaL {wm,Ul)] <TIj . 

Finally we replace each parenthesized term in (|46[) with its corresponding monodromy ma- 
trix, which gives 



M- 



(47) t({u}lAMm) = (-)^*M7^i(i«i,Mi)'iaf ...a^,TM(^M,Mi)'-a 
Letting the monodromy matrices in (j47p be written in the form 

for all 1 < m < Af and contracting on the quantum spaces Vi, . . . , Vm, we recover the result 

62). 

□ 

Lemma 5. The spin-up states | '[\m) and (ffM | are eigenvectors of the diagonal elements 
of the monodromy matrix: 

M M 

(49) A{u, {w}m)\ 1[m) ^X{[u-w,+-i]\ ^m), D{u, {w}m)\ ^m> = !!["- ^^f), 

M M 

(50) (^A/ \A{u,{w}m) = X{[u-Wj+-i]{-^M I, (^M |i?(u, Wa/) = n["~^j]<^^^ I- 

In addition, the spin-down states | ^m) and (JJ-m | are eigenvectors of the diagonal elements 
of the monodromy matrix: 

M M 

[u — Wj + 7] I -IJ-Af), 

M M 

(52) (^A/ \a{u,{w}m) = nt"-^^]^^^^ I' \d{u,{w}m) = n["~^'j+^]<^^^ I- 

2.7. Recovering from the transfer matrix. Let 

(53) t(u, {w}m) = A(m, {w}m) + D{u, {w}m) 

denote the transfer matrix of the XXZ model. The Hamiltonian is recovered via the 
quantum trace identity 

(54) H = [-i]-^\ogt{u) ^, where t{u) = t{u,{w}M) 



Wi = --- = WM- 



Therefore all eigenvectors of t{u,{'w} m) are also eigenvectors of H. The attention turns, 
therefore, to finding vectors € F satisfying 

(55) (a{u, {w}m) + D{u, {w}m)) I*) - T*(u, {wju)]"^) 

for some suitable constants Tq,{u, {w}m). 



2.8. Bethe Ansatz for the eigenvectors. The eigenvectors of t{u,{w}M) are given by 
the Ansatz 



(56) 



1^-) = B{vi,{w}m) ■ ■ ■ B{vN, {w}m)\ 1\m) 



where we assume that N < M, since we annihilate the state | ftA/) when acting with more 
-B-operators than the number of sites in the spin chain. Similarly, we construct eigenvectors 
of t(w, {w}m) in the dual space of states via the Ansatz 



(57) 



{1\m \C{vn, {w}m) ■ • ■ C{vi,{w}m), 



where we again restrict N < M. To ensure that and (|57|) are genuine eigenvectors, 
the variables {ui, . . . ,v^} are required to satisfy the Bethe equations. Using the commu- 
tation relations and the actions ((501 it is possible to show that are 
eigenvectors of t(w, {w}Af) if and only if 



(58) 



M 

n 



\V^ 



7] 



N 

n 



7] 



7] 



for all 1 < I < A'^. Throughout the rest of this work, we will call the expression ([56)) a Bethe 
eigenvector and assume implicitly that the Bethe equations (j58p hold for the parameters 
{v}n- 



3. Domain wall partition function Z 



N 



(^{v}n,{'w}n 



In this section we study Z]s[, the domain wall partition function of the XXZ spin-i chain. 
This quantity acquires its name because it is equal to the partition function of the six- vertex 
model under domain wall boundary conditions [Ij. The calculation of is essential for 
the explicit evaluation of more complicated objects within the XXZ model, such as its scalar 
product. 

3.1. Definition of Zn{{v}n, {w}n)- Let {v}n = {vi, . . . , vn} and {w}n = {wi, . . . , wn} 
be two sets of variables. The domain wall partition function has the algebraic definition 

N 

(59) Zn[{v}n,{w}n) = {U\I[B{u,,{w}n)\1\n), 

1=1 

where the ordering of the _B-operators is irrelevant, since by (1331) they commute. 

3.2. Graphical representation of partition function. Using the graphical conventions 
described in H2.61 the domain wall partition function may be represented as the N x N 
lattice shown in figure [31 



— 1>' 



. f N 



T T T T T 

Figure 3. Domain wall partition function of the six- vertex model. The top 
row of upward pointing arrows corresponds with the state vector | ffAr). The 
bottom row of downward pointing arrows corresponds with the dual state 
vector (Jj-jv I- Each horizontal lattice line corresponds to multiplication by 
a B-operator. 



3.3. Conditions on Zn{{v}n, {w}n)- In [5], Korepin showed that the domain wall parti- 
tion function satisfies a set of four conditions which determine it uniquely. We reproduce 
these facts below, with the fohowing two lemmas. 

Lemma 6. Let us adopt the shorthand = Zn{{v}n,{w}n)- For all TV > 2 we claim 
that 

1. Zn is symmetric in the {w}jv variables. 

2. Zn is a trigonometric polynomial of degree — 1 in the rapidity variable vn- 

3. Setting vn — wn — 7, Zn satisfies the recursion relation 



(60) 



N-1 



wn][wn - Wi- -i]Zn~i, 



I — J. 

where Zn-i is the domain wall partition function on a square lattice of size N — I. 

In addition, we have the supplementary condition 
4. The partition function on the 1x1 lattice is given by Zi — [7] . 
Proof. 

1. We write the domain wall partition function in the form 

(61) Z N (j_v} N , l"^} N )={rN\® {^N \t[{v}n, {w}n) \i[N)®\ Vn) 

with T({v}n, {w}n) given by PO]) . and where we have defined the auxiliary states 



N 



N 



(62) 



N 



t* 



n 



Using the expressions (HTJ, (|i2t for T{{v}n, {'w}n) and contracting on the quantum spaces 
Vi,...,Vn gives 

(63) Zn {{v}nAw}n) ^{Vn \C{wu{v}N)...CiwNAy}N)\rN)- 

The diagrammatic interpretation of the equivalence (j63p is shown in figure [H 



— 1> 

— 1>™1 



T T t T T 



T T T t t 



Figure 4. Equivalent expressions for the domain wall partition function. 
The diagram on the left represents a stacking of the operators B{vi, {w}n)- 
The diagram on the right represents a stacking of the operators C{'Wi, {v}n)- 

By equation ([5^ the C{wi, {v\n) operators all commute, proving that Zn({v}n, {w}n) 
is symmetric in {w}n- 

2. By inserting the set of states Y^^=i \ -D-Jv) (-D-w after the first B-operator appearing 
in (|59p we obtain the expansion 

N N-l 

(64) Zn({v}nAw}n) - Y.^^« \B{vN,{w}N)at\ In){In K \{ B{v,,{w}n)\ 1\n), 



in which all dependence on vn appears in the first factor within the sum. Hence we proceed 
to calculate (JIat \B{vn, {if'}Af)f^| Jlw) for all 1 < n < TV, as shown below: 




Wi Wn Wi Wn W jy 

TTTTT TtTTT 

Figure 5. Peeling away the bottom row of the partition function. The 
diagram on the left represents (JJ-at \B{vn, {w}N)<7n \ ^n), with the internal 
black arrows being summed over all configurations. The diagram on the 
right represents the only surviving configuration. 



The right hand side of figure [5] is simply a product of vertices. Replacing each vertex 
with its corresponding trigonometric weight (see figure 1), we conclude that 

(65) {ii-N \B{vN,{'^}N)o-n\ ^n) = Yl ['"N-Wi][^] J]^ [vn - m + 

l<i<n n<i<N 

Substituting (|65)) into the expansion ((64| gives 

N JV-1 

(66) zn = [7] X! n ^ n [""^ -Wi+^f] 

n=l l<i<n n<i<N i=l 

From this equation we see that every term in Z n {{v} n , {w} n) contains a product of ex- 
actly — 1 trigonometric functions with argument v^. Therefore Z n {{v} n , {w} n) is a 
trigonometric polynomial of degree A'^ — 1 in the variable vn- 

3. We start from the expansion (j66p of the domain wall partition function, and set vn = 
wn — J- This causes all terms in the summation over 1 < n < A'^ to collapse to zero except 
the n = N term, giving 

N-l N-1 

(67) Zn[{v}nAw}n) = [7] it [w^JV - - 7](4jv TT B(wj,{w;}jv)| tTjv)- 



4=1 



We then consider the graphical representation of (JJ-Af |cr^ Ili^i^ ^(wi, {i(;}jv)| ^In), as shown 
below. 



— 1> ► 
— 1> ► 



T 



T T T 



Wi 

T T 



T T 



Figure 6. Peeling away the right-most column of the partition function. 
The left hand side represents the quantity (-IJ-jv Wn Ylfji^ B{vi, {w}n)\ ^In), 
with the internal black arrows being summed over all configurations. The 
diagram on the right contains all surviving configurations. 



The right hand side of figureE]represents the (N—l) x {N— 1) domain wall partition func- 
tion, multiplied by a column of vertices. Replacing these vertices with their trigonometric 
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weights, we find that 

(68) {ii-N Wn n B{v,,{w}n)\1\n) = Y\_\-'"'- " ^n]Zn-i({v}n-i,{'^'}n-i 

i=l 1=1 

Substituting into (p7|) produces the required recursion relation (pO]) . 
4. Specializing the definition (15^1) to the case iV = 1 gives 

(69) ^1(1^1,^1) = (4l \B{vi,Wi)\ ill) =ra,® 1*1 RaAvi,Wl)tl®ia,= [l], 

as required. Alternatively, the lattice representation of Zi is simply the top-right vertex in 
figure m whose weight is equal to [7]. 

□ 

Lemma 7. Let {Ziy}N^fi denote a set of functions Zn{{v}n, {w}n) which satisfy the four 
conditions of the previous lemma. Then Zn = Zat for all iV > 1. In other words, the 
conditions imposed on the domain wall partition function determine it uniquely. 

Proof. From condition 4 on Zi, Z\ we know that Z\ = Z\. Hence we may assume that 
Ztv-i — Zn^i for some N > 2. Using this assumption together with condition 3 on Zn, Zn 
yields 



N-l 

(70) Zn = [7] TT K - wn][wn -Wi- j]Zn-i, 

2—1 

Af-1 

[7] - wn][wn - Wi- "i]Zm-i = Z, 



i=l 



Condition 1 on Zn,Zi^ states that both are symmetric in the variables {wIat. Using this 
fact in the previous equation, we find that 



(71) 



— Zn 



for all 1 < i < A^, 



which proves that Zjv and Zjv are equal at TV distinct values of vn- By condition 2, both 
functions are trigonometric polynomials of degree iV — 1 in u^v, so their equality at N points 
implies Zn = Zn everywhere. This completes the proof by induction. 

□ 

3.4. Determinant expression for Z^. In [9], Izergin found a function which satisfies the 
four conditions of the previous subsection, and is therefore equal to the domain wall partition 
function. We present this formula below. 

Lemma 8. For all > 1 we define 

N 



(72) Zn(^{v}n,{w}n 



n [Vz - Wj + -f][Vi - Wj] 



■det^ W 



l<iJ<N 



l<i<j<N 

det - Wj + l][vk - Wj] 
_ \ / l<iJ<N 

l<i<j"<Af 

The functions {ZyvlArgN satisfy the four conditions of lemma 6. Equivalently, the domain 
wall partition function Zn is equal to the right hand side of (1721) . 

Proof. 
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1. Permuting Wm Wn in dot ( [7] JlfcLi I^fc ~ + t] I^fe ~ ^j] ) swaps two columns 

V ' / l<i,j<N 

of the determinant, which introduces a minus sign into the numerator of (j72p . Similarly, 
permuting Wm ^ Wn in Y\i<i<^j<Ni'^3 ~ introduces a minus sign into the denominator 
of (j72p . These minus signs cancel, leaving Z^r invariant under permutations of the {w}n 
variables. 



2. The numerator of ([72|l is a trigonometric polynomial of degree 2Af — 2 in vn, with zeros 
at the points Vn = Vn for all 1 < rt < TV — 1, since such a substitution would render two 
rows of the determinant equal. The denominator of ([7^ is a trigonometric polynomial of 
degree N — 1 in ujv, with zeros at the same points. Cancelling these common zeros, Zjsi is 
a trigonometric polynomial of degree — 1 in ujv. 

3. Expanding the determinant in Zj^ along its A^"^ row we find 

N 

(73) Zn(^{v}n,{'>^}n^ = Zn-1 f {t;}Ar-i, {wi, . . . , w}, . . . , wn}) 

N-1 N 

k=i 



N-1 N ' 

JI bfc - Wjv] JI [Wj - Wk] 
k=l k=i] 

where Wj denotes the omission of that variable. Setting vn = wn—J in the above expression, 
all terms in the summation over 1 < j < N vanish except the j — N term and we obtain 



(74) Zn = [-f] T\ [vk - wn][wn - Wk - ■y]ZN-i[{v}N-i,h"}N-i 

k—1 

which is the desired recursion relation. 

4. From the definition ([7^ it is clear that Zi = [7]. 

□ 



4. Scalar products Sn{{u}n,{v}N,{'^}Mj 

In this section we define and calculate a sequence of intermediate scalar products Sn, 
which interpolate between the domain wall partition function and the full scalar product. 
The domain wall partition function corresponds to the case n = 0, whereas the full scalar 
product is given by the case n — N . These functions were originally studied in [T3], using 
Drinfel'd twists in the algebraic Bethe Ansatz setting |TS]. Our approach owes much to the 
work of [13 but is somewhat different, since it does not involve constructing the _F-basis. 

4.1. Definition of S'„({m}„, {v}n, {w}m)- Let {u}„ = {ui, . . . , -u„}, {v}n = {wi, . . . , wat}, 
{w}m = {wi, . . . , Wm} be three sets of variables whose cardinalities satisfy < rt < iV and 
1 < < M. We proceed to introduce functions Sn{{u]m{v} n t{w} m) for all < n < A^. 
In the case n = we define 

N 

(75) 'S'o({w}w, {w}m) = (-D-JV/M I W B{vj,{w}m) \ TTm), 

where for conciseness we have suppressed dependence on the set {m}o = 0- As we will soon 
show, up to an overall normalization the scalar product Sq is equal to the domain wall 
partition function Zm- Next, for all 1 < n < A^ — 1 we define 

n N 

(76) Sn{{u}nAv}NAw}M^ = Hn/M I H ^(""^ ' ^ J'^'^) 11 -^('^■J' ' ^J'^'^ ) ' 1^*^^' 

i=l j=l 

12 



where we have adopted the notation N = N — n, which is used frequently hereafter. Finally, 
for n — N we define 



N N 

Sn({u}n,{v}n,{w}m^ = (tTA/ [^^^^(UiJwlM) ]^-B(dj,{w}m)| IIm), 

i=l j=l 



(77) 



which represents the usual scalar product. In all cases ([75)) - ([77)) we shall assume that the 
parameters {v}n obey the Bethe equations (j58p . while the remaining variables {u}n are 
considered free. Accordingly, we name these objects Bethe scalar products. It turns out 
that {>S'n}o<ji<Ar are related by a simple recursion. Hence they provide a convenient way of 
calculating Sn, starting from Izergin's determinant formula (j72l) for Z^- 



4.2. Graphical representation of scalar products. We now provide lattice representa- 
tions of the Bethe scalar products {5'„}o<n<Ar. As in the case of the domain wall partition 
function, the lattices simplify the calculation of these functions. 







— 1> 
— 1> 
— 1> 

— 1> 



N N + 1 

t t t t T T T 







Figure 7. Lattice representation of 5*0. The top row of arrows corresponds 
with the state | ft a/), while the bottom row of arrows corresponds with the 
dual state {ii-N/M I- Each horizontal lattice line represents a S-operator 
B{vj, {w}m)- 



— 1> — ► 




— 1> 

— > ► 



-> ^ 
— 1>" ■* 



1 JV JV + 1 M I if; I 

TTTTTTTTTT 

Figure 8. Lattice representation of Sn- The top row of arrows corresponds 
with the state | ilAj), while the bottom row of arrows corresponds with the 
dual state (JJ-jv/a/ I- '^^'^ highest N horizontal lines represent i3-operators 
B(vj, {w}m)i while the lowest n horizontal lines represent C-operators 
C(mj, {w}m)- 
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-1>1 

— 1> 
— > 

— 1> 



I I I I I I I I I I 

Figure 9. Lattice representation of Sn- The top row of arrows corresponds 
with the state | ffM), while the bottom row of arrows corresponds with the 
dual state (ffM |- The highest N horizontal lines represent i3-operators 
B{vj,{'w}M), while the lowest N horizontal lines represent C-operators 
C{ui, {w}m)- 

4.3. Conditions on S'„({w}„, {v} n, {w}m)- Progressing in the same manner as in WS.'Sl we 
will show that the Bethe scalar products Sni{u}n,{v}N-,{w}M) satisfy a set of properties 
which determine them uniquely. 

Lemma 9. We adopt the usual shorthand Sn — Sn{{u}n, {^'}n, {w}m)- For all 1 < n < 
we claim that 

1. S„ is symmetric in the variables {wj:;j_^^, . . . , wm}- 

2. S'„ is a trigonometric polynomial of degree M — 1 in u„, with zeros occurring at the 
points Un = Wi — 7, for all 1 < j < TV. 

3. Setting u„ — Wj^^^, Sn satisfies the recursion relation 

M 

(78) Sn = TT[w^ I 1 - + Tl'S'n-l, 

where Sn-i denotes the Bethe scalar product S'„_i({u},i_i, {v}n, {w}m)- 
In addition, we have the supplementary condition 

4. 5*0 and Zn are related via the equation 

N M 

(79) -So^l^^liv, {w}m) = ]^ Yl ["'i. - Wj]ZN(^{v}NAw} 

i=l j=N+l 

where we have defined {w}n — {wi, . . . , w^}. 
Proof. 

1. We introduce the auxiliary state vectors 



IN 



N 



N 



(80) (^^|=(g)t:,, (^,'j=(g)4, l^^)=(g)ia., \1ll)=^% 

1=1 2=1 i=l 1=1 

which allow us to write 

(81) Sn - (^^/M l^iVNl® itn \t{{v}n U {u}n, i^} m) \ 1^1) ® \ Vn) ® I ^m), 

where we have defined 

(82) T({t;}ArU{u}„,MM) = 

Tb^{Un, {w}m) ■ • • '7fci(Ml, {w}M)Ta;^{vN,{w}M) ■ --Ta^ivi, {w} m) ■ 
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By application of lemma 4 we thus obtain 

(83) t(^{v}n U {u}„, {w}m) = {-)^'^Ti{wi, {v}n U {u}„) . ..Tm{wm,{v}n U {u}„), 
where for all 1 < i < M we have set 

D{wi,{v}N U {u}n) -B{Wi,{v}N U {u}n 



(84) r,(w„ MjvU {?!}„) = 



-C{Wi,{v}N U {u}„) yl(wi, {wIat U {u}n) 



with {z^Iat U = {vi + 7, . . . , Wat + 7, Ml + 7, . . . , M„ + 7}. Finally, contracting on the 
quantum spaces Vi , . . . , Vm gives 

(85) 



.(j^u}n,{v}N,{w}M 



N 



M 



The diagrammatic interpretation of (|85p is shown in figure [TUl 



— 1> 



— 1> 



>iV+ 1_ 



— > 



TTTTTTTTT 

Figure 10. Alternative graphical representation of S'„. Neglecting an over- 
all minus sign, S„ is equal to this lattice, which is essentially a rotation of 
figure [HI The top row of arrows represents the state | -IJ-^) (g) | tfj^), while 
the bottom row of arrows represents the dual state \ (?) (JJ-Jl |- The 
lowest N horizontal lines represent C-operators, with the remaining lines 
representing _D-operators. 



By equation p9p the D-operators in (|85p commute, proving that Sn is symmetric in the 



variables {i 



'N+l' ■ 



■,Wm}- 



2. Inserting the set of states J2m>N '^ml ^n/m)(^n/m \^rn after the first C-operator ap- 
pearing in (|76p . we obtain the expansion 

(86) 



S'„(^{u}„,{w}jv,{w}m 



N 



^^N/M \C^i^n,{w}Al)<y,n\^N/M)(^ff/M Wni Yi C{u„ {w} m) Yl B{Vj , {w} m)\ i[ 
m>N ^=1 i=l 
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M) 



in which all dependence on u„ appears in the first factor of the sum. We therefore wish to 
calculate (-IJ-jv/m {w}m)o'~ | -U-tv/a/) N < m < M, and do so by identifying it 

with the string of vertices shown below: 



ir ir y A A 

1 JV JV + 1 



= — >' 



JV iV + 1 



TTTTTTTTTT 

Figure 11. Peeling away the bottom row of S'„. The upper diagram rep- 
resents (-IJ-JV/M {w}m)<7~ I -IJ-tv/m)' ■^ith the internal black arrows 
being summed over all configurations. The lower diagram represents the 
only surviving configuration. 



Replacing each vertex in figure [TT] with its corresponding trigonometric weight, we find 
that 



(87) 



JV 



i^N/M {w}M)cr„i\ ^t^/m) = lit"" ^ + 7] J| [u„ - W^][■y] J| [u„ - + j] 

i=l N<i<m m<i<M 

Substituting ([57)1 into the expansion ([55t gives 

JV 

Sn ^ [l]W[Un - Wi + -f] [Un~Wi] [m„ - + 7] 



m>JV N<i<m 
n-1 N 



m<i<M 



1=1 



From the last expression it is apparent that 5„ is a trigonometric polynomial of degree 
AI — 1 in Un- Furthermore, N of the zeros of this polynomial are contained in the factor 



UfLilun - +7]. 

3. Setting u„ — w^j^^ in equation all terms in the sum collapse to zero except the 

term corresponding to m — A'^ + 1 and we obtain 



(89) Sn 



AI 



n-1 



N 



= Ilt^JV+l - ^» +7](4^+i/A./ I n C!iUr,{w}M)Y[B{Vj,{w}M)\ i[ m) , 



i=l 

where we have performed the trivial rearrangements [7] = ["W^jv+i ^ ^jv+i ^ ^Q^^ as 

^^N/M ^ (^JV+i/A/ I produce the first fine of ((89)) . while the second line follows 

directly from the definition of Sn-i- Hence we have proved the recursion relation (j78p . The 
diagrammatic interpretation of this identity is given below. 

16 







— 1> 

— > ► 



1 JV JV + 1 M I i/7 I 

TTTTTTTTTt 

Figure 12. Freezing the last row of the Sn lattice. In general, the vertex 
at the intersection of the u„ and Wj^^-^ lines can be of type 6_(u„, w^^j^) 
or c^{un,Wj^_^^). Setting w„ = Wfj^-^ causes all configurations with a 
b-(um Wf^_^^) vertex to vanish, and we are left with a frozen row of vertices 
as shown. This row of vertices produces the prefactor in ([75)) . whilst the 
remainder of the lattice represents Sn-i- 



4. From the definition of the state vectors ([50)1 we obtain 

(90) So(^{v}n,{w}m ) = {^N/M I {Vn \t{{v}n, {w}m) Wn)®] i^M), 

with T{{v}j\j, {w}m) given by ([3D]). Using lemma 4 and contracting on the quantum spaces 
Vi,...,Vm gives 



JV M 

(91) ^o(Wjv,WAf) =(-)(*^+^)^(mnCK,{^V) n D{w„{v}N)\¥k) 



j=N+l 



Now since | JJ-at) is an eigenvector of the D-operators, as can be seen from equation ([ST]) , we 
have 



JV M JV 

(92) 5o(wjv,Wjv/) = (-)(''+^'"n n K-^»+7](mn^(^-Wiv)m) 



i=l j=JV+l 



or equivalently, substituting Vi — Vi + "f for all 1 < * < into the previous equation. 



(93) 



JV M JV 
1=1 7=JV+1 i=l 



Comparing with the alternative expression (j63p for the domain wall partition function, 
equation completes the proof of (|7^ . A graphical version of this identity is given 
below. 
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— > 



1 N N + 1 M 

TTTTTTTTTt 



Figure 13. Equivalence between and Zjy. The final M — iV columns 
of the Sq lattice must assume the configuration shown. All other configu- 
rations vanish. The block of vertices thus obtained corresponds with the 
prefactor in ([7^ . whilst the remainder of the lattice represents Z^. 



□ 



4.4. Determinant expression for S'„({u}„, {w}jv, {w}m)- 
Lemma 10. Let us define the functions 



(94) f,(n^)^-J:^l[[v,-w + j], 

N M 



N 



(95) g,{u) 



[7] 



TV M 

7] - n [i^fc - " - 7] n [" - wfc] 

k^i 



k=l k^i k=l 

Using these definitions, we construct the N x N matrix 

(fi{wi) ■■■ /i(w^) 9l{Un) 
: : : : 

fN{wi) ■■■ fN{wf}) gN{un) ■■■ 5w(wi) 
Assuming that the parameters {v}n satisfy the Bethe equations (j58l) . we have 

AT M 

n Ilt^^ ~ ^fj] dctMn[{u}n, {v}n, {w}m 

(97) Sr 



n N 

nnt 



n 



n 



n 



l<i<j<n 



l<i<j<N 



l<i<j<N 

The expression (|97p for the intermediate Bethe scalar product Sn originally appeared in 
appendix C of [TB] . 

Proof. Firstly, one must show that Sn is uniquely determined by the set of conditions in 
lemma m This is accomplished using similar arguments to those presented in lemma [71 and 
we shall assume this fact a priori. Hence it will be sufficient to show that the expression 
(l97l) satisfies the list of properties given in lemma [9l 



1. All dependence of the expression (j97|) on the variables {wj^_^^, . . . ,wm} occurs in the 

factor riiLi rijii [^i ~ '^j] in the functions gi{uj) in the determinant. Clearly, these 
terms are invariant under the permutation Wi Wj for all i ^ j. Hence the expression (j97p 
is symmetric in {w^j^-^^, . . . , wm}- 

2. Consider the expression (IMI) for gi{un). Since the variables {w}jv satisfy the Bethe 
equations the numerator of gi(un) vanishes in the limit m„ — >■ Vi. It follows that 
the pole in (j95p is removable, and therefore gi{un) is a trigonometric polynomial of degree 

A^ — 2 in u„. Using this fact, we see that ([W)) is a quotient of trigonometric polynomials 
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in u„. The polynomial in the numerator has degree M + N — 2^ while the polynomial in 
the denominator has degree N — \. We must show that every zero in the denominator is 
cancelled by a zero in the numerator. 

Setting Un = Uj for 1 < j < n—1 causes two columns of the determinant to become equal, 
producing n — I zeros in the numerator which cancel n — 1 of the zeros in the denominator. 
Furthermore since 

hi 



(98) g,{w,) 



Vi - w, 



Y[[vk - Wj + 7] n - + 7] = n ["'j ^ + 7]/»(wj), 



k=l 



k=l 



it follows that by setting w„ = Wj for 1 < j < N , two columns of the determinant are equal 
up to a multiplicative factor, producing N zeros in the numerator which cancel N of the 
zeros in the denominator. This proves that the expression (j97p is a trigonometric polynomial 
of degree M — 1 in m„. 
Finally, since 

[7] ^ ^ 
(99) 9i{wj - 7) = 7— -— r Wlvk - Wj] \Y [wj - Wfe - 7], 

Y"^ + 7J 1^^. 

M N 



k=l 



[Vk ~ Wj + 7] ■ 



k=l k=l 

we see that by setting u„ = Wj ^ "f for all 1 < J < -/V, two columns of the determinant are 
equal up to a multiplicative factor, producing the N zeros which (j97p requires in order to 
satisfy property 2. 

3. Using equation (|98| and the definition of the matrix (|96| . it is clear that 
(100) 

M 



det M„ (^{u}n, {v}n, {w}m^ 

II ["^^+1 -''"fe +7] dctMn-l(^{u}n-l,{'"}N,{w}My 



k=l 



Furthermore, we notice the trivial product identity 

(n N 
i[i[[w,-wj] n n K-^'] 
i=lj=l l<i<j<n l<i<j<N 



n-1 N+1 

II II ^ ^3 



n n 

1=1 j=l l<i<j<n-l l<i<j<N+l 

Combining the results (llOOp and (|101[) . we find that the expression (P7|) satisfies the recursion 
relation 

4. Taking the n = case of (|97|) yields 



N M 



(102) 



II nt"^ -WjldetA^o 



l<i<j<N 

with the matrix AIq given by 

/ j:;j^=^Uk^li^k-Wi+-f] ••• 



(103) Mo 



[7] 



rifc^^iK - i"iv +7] \ 



JiL 



■ri/c^jvK -^i'w +7] / 



Comparing (I102p with the determinant expression ((7^ for the domain wall partition func- 

rM 

□ 



tion, we see that Sq = YiiLi IltlAf+ibi ~ Wj]ZN, as required. 
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4.5. Evaluation of Sn{{u}n, {v}]\f, {w}m)- Let us now consider the n — N case of equa- 
tion (|97p in more detail. For purely aesthetic purposes, we simultaneously reverse the order 
of the columns in the matrix M n and the order of the variables in the Vandermonde 
ni<i<j<Af [""i ~ ""il- We also take the transpose of the matrix Mn- The formula (j97p is 
invariant under these transformations, and we obtain 

N M 
1=1 ]=1 



(104) Sn[{u}nAv}nAw}m) ^ 



n K 

l<i<j<N 



det 



/ ^ 



M 



7] lit 

fe=i 



Wk\ 



N 



M 



- 7] W ["i - Wfe + 7] 



The expression (|104l) for the Bethe scalar product was discovered by Slavnov in [2^. The 
original proof required a recursion relation between scalar products of dimension N and 
— 1, which can be found in section 3, chapter IX of jT|. 

5. Trigonometric Felderhof model 

In this section we devote our attention to the trigonometric Felderhof model. Most of the 
material that we present originally appeared in |12| . 

5.1. Hamiltonian H. The Hamiltonian of the XX spin-^ chain in an external magnetic 
field h is given by 

M 

(105) ^ = 2 ^ («+i + + 2/i< 

m— 1 

assuming the usual periodicity of the Pauli matrices. Similarly to the anisotropy constant 
in the XXZ model, we will find it convenient to parametrize h = cosh(2p), where p is a 
variable to be introduced below. 

5.2. _R- matrix and Yang-Baxter equation. At the free fermion point 7 = -y, the R- 
matrix (j20p may be generalized to include extra variables, in such a way that the Yang-Baxter 
equation remains satisfied. This leads to the trigonometric limit of the model introduced by 
Felderhof in [10] , and accordingly we call it the trigonometric Felderhof model. This model 
was also studied in |llj . as the first in a hierarchy of vertex models with increasing spin. 
The i?-matrix for the trigonometric Felderhof model is given by 



(106) Rab{u,p,v,q) ^ 



( a+{u,p,v,q) 



V 





b+{u,p,v,q) 
c^{u,p,v,q) 






c+{u,p,v,q) 
b-{u,p,v,q) 




a-{u,p,v,q) j ^ 



where we have defined the functions 



(107) a±(M,p,w,(7) = [±(u - t;) +P + g], 

(108) h±{u,p,v,q) = \u-v±{q-p% 

(109) c±(u,p,«,(z) = [2p]^[2<z]^ 

with \u\ = sinhu as usualQ The i?-matrix is an element of End(VQ ® Vf,), and the variables 
u, V are rapidities associated to the respective vector spaces Va,Vb. The new features of this 
i?-matrix are the variables p,q. These are called external fields, and are associated to the 
respective vector spaces Va,Vf,- We recover the free fermion point of the six- vertex model 
by setting p ^ q 



7Tl 

4 ■ 



"'^The parametrization of 11' is recovered by multiplying all weights by e" "+P+9 and setting e 
= B. 



2p 
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The entries of the i?-matrix (|106l) admit the same graphical representation as those of the 
i?- matrix (|20p. The only difference is that each vertex line now accommodates a rapidity 
variable and an external field. Hence we identify the functions (|107p - (jl09p with the vertices 
shown below. 



_I>«,P ^ 



a+{u,p,v,q) 



-|>«,p 



b+{u,p,v,q) 



c+{u,p,v,q) 



^ ^ ^ 

a-{u,p,v,q) b^{u,p,v,q) C-{u,p,v,q) 

Figure 14. Six vertices of the trigonometric Feldcrhof model. 



With the following result we assert that the Yang-Baxter equation continues to hold, even 
in the presence of the external fields. 



Lemma 11. The i?-niatrix (jl06p obeys the Yang-Baxter equation 
(110) 

Rab{u,p, V, q)Rac{u,p, w, r)Rbc{v, q, w, r) = Rbc{v, q, w, r)Rac{u,p, w, r)Rab{u,p, v, q). 
This is an identity in End(Va ®Vb® Vc), true for all u, w, w and p, q, r. 










3b 




3a 



^ «,p 





3a 

k 


ia 






k 

ui, r 




T T 

Figure 15. Yang-Baxter equation for the trigonometric Felderhof model. 
Each index ia,'i-b,ic,jaijbijc represents a black arrow that is fixed on both 
sides of the equation, while the indices ka^kb, kc are summed over all possible 
arrow configurations. 



5.3. Monodromy matrix and intertwining equation. The monodromy matrix is an 
ordered product of _R-matrices, given by 

(111) Ta{u,p,{w,r}M) = Ral{u,p,Wi,ri) . . . RaM{u,p,WM,rM), 

with the multiplication taken in the space End(Va). We write the contribution from the 
space End(Va) explicitly, by defining 

ni9^ T Un \-( ^("'P'{^'^}m) B{u,p,{w,r}M) \ 

\ C[u,p,{w,r}M) D{u,p,{w,r}M) J J 

where the matrix entries are all operators acting hi Vi® ■ ■ ■ ® Vm ■ For a graphical repre- 
sentation of these operators, we refer the reader to figure [21 The correspondence is exactly 
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the same, except that the rapidity u is now accompanied by the external field p, and each 
Wj by an rj . Due to the Yang-Baxter equation (jllOp , we obtain the intertwining equation 



(113) Rab{u,p,v,q)Ta{u,p,{w,r}M)Tb{v,q,{w,r}M) = 
Tb{v, q, {w, r}M)Ta{u,p, {w, r}M)Rab{u,p, v, q). 

As usual, this leads to sixteen commutation relations amongst the entries of the monodromy 
matrix (I112p . Of these commutation relations, two have particular significance in our later 
calculations. They are given by 

(114) [u - V + p + q]B{u,p,{w,r}M)B{v,q,{w,r}M) = 
[v-u + p + q]B{v, q, {w, r}M)B{u,p, {w, ^Im), 

(115) [v-u + p + q]C{u,p, {w, r}M)C{v, q, {w, r}A/) = 
[u-v+p + q]C{v, q, {w, r}M)C{u,p, {w, r}M)- 

5.4. Recovering H from the transfer matrix. Let 

(116) t{u,p,{w,r}M) = A{u,p,{w,r}M) + D{u,p, {w,r}M) 

denote the transfer matrix of the trigonometric Felderhof model. The Hamiltonian (|105p is 
recovered via the formula 



d 

(117) H = [2p]—\ogt{u,p) , where t{u,p) ^ t{u,p,{w,r}M) 

OU u—O 



Wi — '-' — WM—O' 

ri = ---=rM=P 



Therefore all eigenvectors of t(M,p, {w, tIm) are also eigenvectors of H. Once again, the 
problem of diagonalizing H is reduced to finding vectors j^") € ^ satisfying 

(118) {a{u,p, {w, r}M) + D{u,p, {w, tIm)) |*) = T*(u,p, {w, r}A/)|*) 
for some constants t*(m,p, {w,r}M)- 

5.5. Bethe Ansatz for the eigenvectors. We construct the Bethe eigenvectors of the 
trigonometric Felderhof model in an analogous fashion to in i i2.8l The eigenvectors 14") S 
and (\E'| G V* of t{u,p, {w,r}M) are given by the Ansatze 

(119) l^-) = B{vi,qi,{w,r}M)...B{vN,qN,{w,r}M)\ IIm), 

(120) (^l = (ftM \C{vN,qN,{w,r}M)---C{vi,qi,{w,r}M), 

with < M in both cases. To show that these are genuine eigenvectors, it is necessary to 
act on them with t{u,p, {w,r}M) using the commutation relations contained in (|113|) . We 
find that (|119p and (|120p are eigenstates of the transfer matrix if and only if 

M M 

(121) (-)^ Y[ [vj + qj -Wk+rk] + Y[ [v-i + qj - Wk - rk] ^ 

k=l k=l 

for all 1 < j < A^. The identities p2ip constitute the Bethe equations for the trigono- 
metric Felderhof model and in contrast to those of the XXZ chain (|58p , the dependence on 
{«!, . . . ,vn} is decoupled. This is a consequence of the free fermionic nature of the model 
being studied. Hence we may interpret {vi + qi, . . . ,vn + qN} as being roots, in x, of the 
single equation 

M M 

(122) i-f l[[x- Wk +rk] + l[[x-Wk- rk] = 0. 

k=l fe=l 
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N({v,q}N,{w,r}N 



6. Domain wall partition function Z 

6.1. Definition of Zn{{v, q}N, {w, r} n). The domain wall partition function of the trigono- 
metric Felderhof model has the algebraic definition 



(123) 



N 

ZN{{v,q}N,{w,r}N'j = {ii-N \^ B{vj,qj,{w,r}N)\ 1\n) 



This naturally extends the definition of the domain wall partition function (j59p to a model 
containing external fields. Notice that we must define an ordering of the i3-operators in 
(|123p . since by (I114p they do not commute. 

Similarly to before, we represent the domain wall partition function by an x TV lattice, 
as shown in figure 1161 



— > 
— > 

. <3JV - 



T T T T T 

Figure 16. Domain wall partition function of the trigonometric Felderhof 
model. The top row of arrows corresponds with the state vector | ft at). The 
bottom row of arrows corresponds with the dual state vector (JJ-at |. Each 
horizontal lattice line corresponds to multiplication by a B{vj,qj^ {w^r}N) 
operator. Notice that the ordering of these lattice lines respects the ordering 
of S-operators defined in (I123p . 



6.2. Conditions on Z7v({w, 9}jv, {w, r}jv)- We progress towards calculating the domain 
wall partition function (I123|) . The procedure begins with the following result from [T^], 
which establishes a set of Korepin-type conditions on Zn{{v, q} jv, {w, r}j^). 

Lemma 12. We adopt the shorthand Z^ = Zpf{{v, q}N, {w, r}i\f). For all > 2 we claim 
that 

1. Zn is a trigonometric polynomial of degree A'^ — 1 in the rapidity variable vn- 

2. Zn has zeros at the points vn ^ Vj + qj + qN, for all 1 < J < ~ 1- 

3. Setting vn — wn + Qn + fN, satisfies the recursion relation 
(124) 



N-1 



= [2(jAr] 2 [2rN] 2 Y[ [^JV - Wj + Tj + tn] [vj - wn + qj - rNjZ^-i, 



where Zj^^i is the domain wall partition function on a square lattice of size iV — 1. 

In addition, we have the supplementary condition 
4. The partition function on the 1x1 lattice is given by Zi — [2qi]^ [2ri]^ . 
Proof. 
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1. By inserting the set of states J2n=i ^n}{^n |f„ after the first B-operator appearing 
in (|123p . we obtain the expansion 

N 

(125) ZN({v,q}N,{w,r}N'j = ^(^TV |-B(w7v, gjv, {w, r}7v)cr+| iL^v) 

N-1 

X {^N B{vj,qj,{w,r}N)\1\N), 

in which aU dependence on ujv appears in the first factor within the sum. Hence we shah 
calculate {]}-n \B{vN,qN, {w,r}N)a^\ ^at) for all 1 < n < TV, as shown below: 



TTTTT TTTTT 

Figure 17. Peeling away the bottom row of the trigonometric Felderhof 
partition function. The diagram on the left hand side represents the quan- 
tity (JJ-Af \B{vN,q]\f,{w,r}pf)a^\ JJ-Af); with the internal black arrows being 
summed over all configurations. The diagram on the right represents the 
only surviving configuration. 



The right hand side of figure [T7] represents a product of vertices. Replacing each vertex 
with its corresponding trigonometric weight (see figure 14), we have 

(126) {i^N \B{vN,qN,{w,r}N)<J+\M = 

[2(7 at] ^ [2r„] i - Wj + Vj - qN] Y\ I'^'j - ""N + qN + r-j]. 

l<j<n n<j<N 

Substituting (|126p into the expansion (|125p gives 

N 

(127) [2^^] ^ P*^"] ^ n [""^ ~ ^1 + ''j ~ n ["'j -""N + qN + rj] 

71=1 l<j<n n<j<N 

N-1 

X i^N krVfJ B{Vj,qj,{w,r}N)\1\N)- 

From (jl27p we see that every term in Z]y{{v,q}N, {w,r}N) contains a product of — 1 
trigonometric functions with argument vn. Thus Z]\j{{v,q}N,{w,r}N) is a trigonometric 
polynomial of degree A'^ — 1 in the variable . 

2. We multiply the partition function (jl23p by 11^^^ [^n — Vj + qj + qN] and repeatedly use 
the commutation relation 

(128) [vN -Vj + qj + qN]B{vN,qN,{w,r}N)B{vj,qj,{w,r}N) = 
[vj "VN + qj + qN]B{vj,qj,{w, r}N)B{vN,qN, {w, r}Ar), 

which is a rewriting of (|114p . to change the order of the -B-operators. We obtain 

N-l 

(129) Y\_ [""N - Vj + q-j + qN]ZN [{v, q}N, {w, r}N^ = 
i=i 

Af-l AT-l 

W [vj - VN + qj +qN]{ii-N l-ff B{vj,qj,{w,r}N)B{vN,qN,{w,r}N)\ 1\n)- 
Graphically, we depict (I129p with the following diagrams: 
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—> qi 



vn 

qN~ 



ri 



t t t 



— 1> qi 



qN ri 



T T T T T 



Figure 18. Reordering the lattice lines of the trigonometric Felderhof par- 
tition function. The diagram on the left is the domain wall partition func- 
tion multiplied by the string of vertices Y\^^^ a-{vj,qj,VN ,(In), and it 
corresponds with the left hand side of (|129p . Each vertex can be threaded 
through the lattice using the Yang-Baxter equation, which ultimately pro- 
duces the diagram on the right. This diagram represents the domain wall 
partition function with its iV"^ row transferred to the top of the lattice, 
multiplied by the string of vertices J^^"^^ a+(wj, , ujv, ^at). Clearly, this 
corresponds with the right hand side of (I129p . 



The right hand side of (jl29p is a trigonometric polynomial of degree 2 A'^ — 2 in wjv , with 
zeros at the points vn = Vj +qj + qN for all 1 < j < N —1. Therefore the partition function 
Zjv({w, g}Ar, {w, r}Ar) must have zeros at the same points. 

3. We start from the expansion (jl27p of the domain wall partition function, and set vn ~ 
wn + + tn ■ This causes all terms in the summation over 1 < rt < to collapse to zero 
except the n = N term, and we obtain 



(130) ZN(^{v,q}N,{'w,r}N^ 



VN^WN+qN+rN 

N-1 N-1 

[2gAr]5[2rAr]^ J| [wn - Wj + + rN]{i)-N Wn^^ B{vj , qj , {w,r} n)\ ilN)- 
j=i i=i 
We further simplify the previous expression by using the diagrammatic representation of 
Vj,qj,{w,r}N)\ 1\n), shown below: 



TTTTT TTTTT 

Figure 19. Peeling the right-most column of the trigonometric Felderhof 
partition function. The diagram on the left hand side represents the quan- 
tity (JJ-AT Ic^^fl^i^ B{vj, qj, {w, r}N)\ ft at), with the internal black arrows 
being summed over all configurations. The diagram on the right contains 
all surviving configurations. 



The right hand side of figure [TO] represents the (A'' — 1) x {N — 1) partition function, 
multiplied by a column of vertices. Replacing these vertices with their trigonometric weights, 
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we have 

N-1 N-1 

(131) (-ILjv IcTjv^ S(uj,qj,{u;, r}Ar)| tTAf) = ^[[""3 - + Qj ~ rN]ZN-i- 

Substituting (|13ip into p30l) we recover the required recursion relation (|124p . 
4. SpeciaUzing the definition (I123P to the case iV = 1 gives 

(132) Zi{vi,qi,wi,ri) - (^1 |S(«i, gi, {u;, r}i)| ^1), 

^tli i?aii(vi,gi,wi,ri)ti ®iai= [2gi]^[2ri]^ 

as required. Alternatively, the 1x1 partition function is the top-right vertex in figure 1141 
whose weight is equal to [2(7i] 2 [2ri] 2 . 

□ 

6.3. Factorized expression for Zn{{v, qJn, {w, r}]^). The conditions 1-4 are strong con- 
straints. Not only do they specify Zj\j{{v,q}js[,{w,r}N) uniquely, they lead to its direct 
evaluation, as we demonstrate below. 

Lemma 13. The domain wall partition function has the factorized expression 
(133) 

N 

ZN(^{v,q}N,{'^i'>^}N^ = ]^[2gj]^[2rj]3 [y^ ^ y^, ^ ^ q^,^[^^, _ ^ + Vk]. 

j=l l<j<k<N 

The result (|133p was first obtained in jlG^ using a complicated recursion relation. A more 
straightforward proof, based on solving the conditions 1-4, subsequently appeared in |12] . 
It is the latter proof which we present below. 

Proof. From condition 1 and 2 on Zn{{v, q}N, {w, J'} at) we know that it must have the form 

N-1 

(134) Z7v(^{w, g}jv,{«^,r}jv) =c(^{w}jv-i,{g}A',{w, r}Ar) K - + + ^at] , 

i=i 

where C does not depend on w^v, but depends on all other variables. Evaluating (|134p at 
vn = wn + qN + tj^ and comparing with condition 3 on Zjv, we obtain 

(135) 

N-1 

Zn 

ti ht —m -\-n i\r -\-r \ / 

N-1 

= [2qN] 2 [2rN] ^ J| [wjv - Wj + rj + rjv] [vj - wn + q^ - rN]ZN- 
from which we extract the equation 

N-1 

(136) C = [2qN]^2rN]^ [wat - -f -t- r^jZyv-i (^{w, glw-i, {w, rlAr-i^ 

Substituting this expression for C into (jl34p . wc obtain the recurrence 

(137) Z^v({w,g}jv,{w,r}w) = [2qN]H2rN]ix 

N-1 

Y[ [vj -VN + qj +qN][wN -Wj +rj + rN]ZN-i(^{v, q}N-i, {w^'>^}n-i^ , 



= C({t;}Ar_i,{g}Ar,{w,r}Ar) TT [Vj - Wn + q^ - Tat], 



whose basis is given by condition 4. This recurrence is trivially solved to produce the formula 
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n (^{U,p}n, {v, q}N, {W, r} M 



7. Scalar products S, 

7.1. Definition of Sn{{u,p}n,{v,q}N,{w,r}M)- Define three sets of rapidity variables 
{u}n = {ui, ■ . ■ ,Un},{v}N = {vi, . ■ . , vn} , {w} M = {wi, . . . ,wm} and their correspond- 
ing sets of external fields {p}n = {pi, . . ■ ,p„}, {q}N = {qi, . . . ,qN},W}M = {ri,...,rM}. 
The cardinalities of these sets are assumed to satisfy < n < and 1 < iV < M. For n = 
we define 

N 

(138) 5'o({f,9}7V,{w,r}M) = {ii-N/M \^ B{vk,qk,{w,r}M)\ i\M)- 

fe=i 

Similarly to in §4 we will find that 5*0 is equal to the trigonometric Felderhof partition 
function Zn, up to an overall normalization. Next, for all 1 < n < iV — 1 we define 



(139) 

n N 
'5'n({'",_p}„,{-y,g}Ar,{w,r}A/) = i^N/M \^ ("^j ^ Pj A^, r} B{vk, qkA^, r} m)\ 1\ M ) 



J=l fc=l 

with N = N — n. Finally, in the case n = TV we fix 
(140) 

TV N 

SN{{u,p}N,{'",q}N,{'>^ir}M'j = iilM \^C'{uj,pj,{w,r}M)^ B{vk,qkAw,r}M)\1\M)■ 
J=l k=i 

The scalar products (|138p - (|140p are the trigonometric Felderhof analogues of those defined 
in H4.ll They have identical graphical representations to those described in ^4.21 except that 
every rapidity variable is now accompanied by an appropriate external field. In the following 
subsection we give a set of conditions on these scalar products, using similar techniques to 
those developed earlier. 

7.2. Conditions on ^^^({u,^},!, {w, g}Ar, {w, tIm)- 
Lemma 14. For all 1 < n < iV we claim that 

1. Sn is invariant under the simultaneous permutation of variables {wj,rj} {wk,rk} for 
all j,/ce {iV+l,...,M}. 

2. Sn is a trigonometric polynomial of degree M — 1 in u„, with zeros occurring at the 
points Un = Pn + Wj + fj, for all 1 < j < N . 

3. Setting u„ + Pn — ^tv+i + ''tv+i' satisfies the recursion relation 

(141) Sn = [2p„]^[2rj^ ]5 Yl K - -f^jv+i + ''j - ^W+l + 2Pn] 

N+l<j<M 

where we have abbreviated Sn-i = 'S'„_i({u,p}„_i, {u, g}jv, {w, tIm)- 
In addition, we have the supplementary condition 

4. 5*0 and Zj^ are related via the equation 

N M 

(142) So(^{v,q}NAw,r}M^ =Y]_ H K - ""'fc + 9j - '^fcl^w glw, {w, -^1^^ 

j=l k=N+l 

Proof. The proof of properties 1-4 is analogous to the proof of lemma 9. There, we presented 
an algebraic proof of the properties. Here, we outline a less technical graphical proof. 

1. For any + 1 < j < M, multiplying S'„({m,p}„, {w, qIat, {w, r}Af) by the weight 
a^(wj,rj ,Wj^i,rj^i) is equivalent to attaching an a+ vertex at the base of the lattice, 
as shown in figure [201 
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v,q 



u,p 



N 



► 
► 



—0 



N N + 1 



TTTTTTTTTT 

Figure 20. Attaching an a^{'Wj,rj,Wj-i,rj-i) vertex to the Sn lattice. 
The points marked with grey arrows are considered to be summed over all 
arrow configurations, but the only non-zero configuration is the one shown. 



The attached vertex can be translated vertically through the lattice using the graphical 
version of the Yang-Baxter equation, as given by figure 1151 It ultimately emerges from the 
top of the lattice, still as an a+{wj,rj,Wj-i,rj-i) vertex, and the (j — 1)**^ and j'^ lattice 
columns are swapped in the process. The result of this procedure is shown in figure [5TJ 



v,q 



T T T 



M w, r 



Figure 21. Extracting the a+{wj^rj^Wj^i^rj^i) vertex from the Sn lat- 
tice. Once again, the grey arrows indicate the only surviving configuration 
in the summation at those points. 



Cancelling the common factor a+ [wj , Vj ,Wj-i,rj^i) from figures HD] and [5T1 we see that Sn 
is invariant under swapping the (j — 1)"^ and j*"^ lattice columns, for all 7V-f 1 < j < M. An 
arbitrary permutation of the lattice columns is just a composition of such swaps. Therefore 
Sn is invariant under permuting its j"^ and fc"^ columns, for all -I- 1 < j, fc < M. 

2. Consider the graphical representation of the scalar product Sm as given below: 
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v,q 



— 1> 



— > 
-t>' 



JV JV + 1 
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Figure 22. Lattice representation of Sn, with frozen vertices included. 
The grey arrows indicate points which are summed over all configurations, 
but whose only non-zero configuration is the one shown. 



We examine the final row of this lattice, through which the variable m„ flows. Every 
non-zero configuration of this row contains a c-{un,PmWj,rj) vertex, which by (jl09p does 
not depend on m„, and A/ — 1 other vertices which are trigonometric polynomials of degree 
1 in Un- ft follows that Sn is a trigonometric polynomial of degree M — 1 in u„. 

Furthermore, all surviving configurations of the final row contain the N vertices as shown 
in figure mi Consequentially, Sn contains the factor 



JV 



JV 



(143) 



Y[a-{un,Pn,Wj,rj) = Y[h 



which gives rise to zeros at u„ 



Wi 



for all 1 < j < iV. 



3. Consider the vertex at the intersection of the u„ and Wj^^^ lines in figure [22l In any given 
lattice configuration, this can be of type 5_(u„,p„,?ii^ 



Pn 



• j^j^^ cancels all terms containing 6_(u„,p„, 



Setting Un 

freezes the entire final row of the lattice to the configuration below: 



li^jV+l) (^m '^AJ -1-1 ' ''jV-t-l)' 



V,q 



—V 



— > 



JV 



. 1 



— 1> ■* 



1 JV JV + 1 M 

TTTTTTTTTT 

Figure 23. Freezing the entire last row of the Sn lattice. The last row of 
vertices produces the prefactor in (I14ip . while the remainder of the lattice 
represents Sn-i- 



From the diagram we see that setting u„ + Pn — w^^v+i ''tv+i reduces Sn to Sn-i, up 
to a multiphcative factor. This factor is evaluated by matching each vertex in the final row 
with its trigonometric weight, giving 
(144) 



l<j<N+l 

N+l<j<M 

Using the explicit formulae (I107P and (|109l) for the functions appearing in (I144p . we obtain 
the required recursion relation (|14ip . 

4. The scalar product Sq is represented by the lattice below: 



v,q 



— > 
— 1> 
— > 

— 1> 



N JV + 1 



ttttttTTTt 

Figure 24. Frozen vertices within 5*0. The final M — N columns of vertices 
produce the prefactor in (|142l) , while the remainder of the lattice represents 
Zn- 



The vertices in the last M — N columns of a lattice configuration must be of the form 
b-{vj,qj,Wk,rk), or else the configuration vanishes. Peeling away this block of frozen ver- 
tices, we find that 5*0 is equal to Zn up to the overall factor Y[f=i IlfeLjv+i ^- (^i' Ij^'^k, ^k)- 

□ 

7.3. Factorized expression for Sn{{u,p}n,{v,q}N,{w,r}M)- 

Lemma 15. Let vj + qj satisfy the equation (jl21l) for all 1 < j < A^. In the presence of 
this constraint, the scalar product Sn has the factorized expression 

(145) 

n N N 

Sn(^{u,p}n,{v,q}N,{w,r}M) = l[[2pj]^ Y[[2q,]^ Y[[2r,]i X 

j=i i=i i=i 

Y[ [uk - Uj + pj + pk] Yl + qj + qk] Y\. ['^'^ ~ + + ^ 

l<j<fc<ri l<j<fe<7V l<j<fe<iV 

n N N M 

H n [wfe - Uj + Pj +rk]Y\_ H -Wk + qj - rk] X 

J = l J = l k=N + l 

n N 

nn 



1 



M 



M 



,X{[{-rX{[uj+Pj-Wk+rk] + \{[ 

' i=l \ k=l k=l 



+Pi -Wk- Tk] 



1 fc=i ^"j' + -'"k-qk 

Proof. We begin by stating that the conditions 1-4 are uniquely determining^ Hence we 
need only verify that (I145p satisfies properties 1-4. 



This is proved along very similar lines to lemma 7. 
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1. By studying (|145l) we see that Sn has dependence on {wj^^^, . . . , wm} and {r^^^^' • • ■ ' ''J^^i 
only through the terms 

N M n / M M \ 

n n ~ ■^fe + ^3 ^ ^k]-, W {-)^ W [uj + Pj - Wk + rk] + Y\_ ["i + Pj -Wk~ Tk] 

J = l k=N+l J = l ^ '^■=1 / 

which are clearly both invariant under the permutation {wj,rj} o {wk,rk} for all j,k G 
{7V+1,...,A/}. 

2. Since vj + qj is a root of the equation (I12ip for all 1 < j < iV, it follows that 

n N ^ n / M M \ 

n n [u,+p,-v,-q,] n (-)" nt". +^.- - + '^^i + n^- - - -^-i 

j = l k=l ^ 3 ' fJ K ':i'-J \ / 

is a trigonometric polynomial in u„ of degree M — N. The remaining terms in (jl45p comprise 
a trigonometric polynomial of degree — 1 in u„. Therefore the entire expression (|145p 
is a trigonometric polynomial of degree A/ — 1 in w„. In addition, the required factor 

[^j — + Pn + ''j] is present in (jl45p . 

3. The recursion relation (|14ip is proved by setting u„ + p„ = w^j^-^ + r^_^^ in p45p 
and rearranging the factors in the resulting equation. Since this procedure is tedious but 
elementary in nature, we shall omit the details. 

4. Setting n = in (|145p gives 

N M 

(146) So(^{v,q}N,{w,r}M^ =Yl H ["j - '^k + Qj - rk]x 

j=l k=N+l 

N 

Y\_[2qj]^[2rj]^ Y[ [vj - Vk + qj + qk][wk - Wj + rj + rk]. 

j=l l<j<k<N 

Comparing equation p46p with the factorized expression p33p for the domain wall partition 
function, we verify (|142p . 

□ 

7.4. Evaluation of SN{{u,p}N,{v,q}N,{w,r}M)- For completeness, we write the n — N 
case of equation (|145p explicitly. We have 

(147) 



TV 

^ n PPj] ^ 2 Y[ [uk-u-j + Pj + pk] [vj -Vk + qj + qk] x 

j=l i<j<k<N 
N M 

UU^'^.-Wk+q.-rk] ^ ^ 

n ( n ["j' + -'^k+rk] + Y[ K + Pj - '^k -rk] \ ■ 

nr , 1 3 = 1 \ k=l k=l ) 

m + Pj -'"k- qk\ 

j,k=l 

Let us consider the zeros of this expression in the variable +pn- Studying the first line 
of the right hand side, we see that TV — 1 of the zeros are of the same type as in the domain 
wall partition function. The remaining zeros, in the second line, coincide with the M ~ N 
roots of the Bethe equation (jl22p which are different from {vi + qi, . . . ,vn + q^}- 
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